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Abstract 

Sequences based on the Distant Dipolar Field (DDF) liave shown great promise for novel 



^' spectroscopy and imaging. Unless spatial variation in the longitudinal magnetization, M^{s), 

(— I . is eliminated by relaxation, diffusion, or spoiling techniques by the end of a single repeti- 

^ Ph | tion, unexpected results can be obtained due to spatial harmonics in the steady state Mf^{s) 

profile. This is true even in a homogeneous single-component sample. We have developed 
an analytical expression for the Mf^{s) profile that occurs in DDF sequences when smear- 
ing by diffusion is negligible in the TR period. The expression has been verified by directly 
Tj- ' imaging the Mf^{s) profile after establishing the steady state. 
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1 Introduction 



NMR and MRI sequences utilizing the Distant Dipolar Field (DDF) have the rela- 
tively unique property of preparing, utilizing, and leaving spatially-modulated lon- 
gitudinal magnetization, My(s) , where s is in the direction of an applied gradient. In 
fact this is fundamental to producing the novel "multiple spin echo"[l, 2] or "non- 
linear stimulated echo" psl] of the classical picture and making the "intermolecular 
multiple quantum coherence (iMQC)" ||4J] observable in the quantum picture. 
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Figure 1. Pulse Sequence. All RF pulses shown as hard for simplicity are actually Sinc3. a 
and P are the same phase. 

Existing analytical signal equations for DDF/iMQC sequences depend on M^{s) 
being sinusoidal during the signal build period] 5, 6]. Experiments that probe sam- 
ple structure also require a well-defined "correlation distance" which is defined as 
the repetition distance of My{s) I2ll8i,|9|]. If the repetition time TR of the DDF se- 
quence is such that full relaxation is not allowed to proceed TR < 5Ti,ot diffusion 
does not average out the modulation, spatially-modulated longitudinal magnetiza- 
tion will be left at the end of one iteration of the sequence. The next repetition of 
the sequence will begin to establish "harmonics" in what is desired to be a purely 
sinusoidal modulation pattern. Eventually a steady state is established, potentially 
departing significantly from a pure sinusoid. 



2 Experimental Methods 



In order to study the behavior of the steady state M^^{s) profile we have imple- 
mented a looped DDF preparation subsequence followed by a standard multiple- 
phase encode imaging sub-sequence. (Figure 1.) The a pulse excites the system, 
the gradient Gq twists the transverse magnetization into a helix. (3 rotates one com- 
ponent of the helix back into the longitudinal direction. For simplicity we have 
omitted the 180° pulses used to create a spin echo during TM and/or TB sometimes 
present in DDF sequences. Also, we are only interested in Mz{s) in this experi- 
ment, not the actual DDE-generated transverse signal. Looping the "preparation" 
sub-sequence thus creates the periodic M^{s) profile, spoils remaining transverse 
magnetization, and establishes Mf^{s). The £ pulse converts Mf^{s) into transverse 
magnetization, allowing it to be imaged via the subsequent spin echo "image" 
sub-sequence. M^^{s) must be re-established by the "preparation" sub-sequence 
for each phase encode. After a suitably long full relaxation delay "relax," the se- 
quence is repeated to acquire the next k-space line. This is clearly a slow acquisition 




Figure 2. Theoretical values of M^{s). M^^{s)is shown dashed 
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envelope, M^ {s) is shown as a heavy line, M^ ' {s) as a normal line, 
a = p = 90°, r/? = 2s, TM = 0ms, TB = \00ms, Ti = I As 

method because many TR periods are required to reach steady state in the prepa- 
ration before each k-space line is acquired. The sequence is intended as a tool to 
directly image the M^^{s) profile, verifying the Mf^{s) that would occur in a steady 
state DDF sequence, not as a new imaging modality. 



3 Theory 



The effect of the "preparation" pulse sequence was first determined for a single 
iteration. The progress along the sequence is denoted by the the superscript. 



Starting with fully relaxed equilibrium magnetization before the a pulse: 

after the a pulse, the mix delay TM and the |3 pulse we have: 

M^{s) = [A^cos{qs) +B^]M^'' + C^Mo 



TM 



(1) 



(2) 



aP = -sin{a) e ^2 sin{^) 
5P = cos{a)e ^1 cos{^) 

CP = (l-e ^1 )co5(p) 

The parameter q = ^, where X is the helix pitch resulting from the applied gra- 
dient. Diffusion has been assumed to be negligible at the scale of X. Note that T2 
is used in A rather than Tj when Gq is larger than background inhomogeneity and 
susceptibility gradients. 



After the build delay TB we have: 



Ml^{s) = [A^^cos{qs) +B^^]M^''{s)+C^^Mq (3) 

_TM_ _TB 

A^ ^ = -sin{a) e ^2 sin(^) e ^1 

_TB TB 

B^^ = cos{a)e ^ico5(p)e ^1 
C^« = [(l-e ^i)co5((3)-l]e ^1 + 1 

At the start of the next repetition, after a TR period inclusive of TM and TB we 
have 

mJ^{s) = [A'^'^cosiqs) +5^^]Mf«(5) +C^'^Mo (4) 

_7;m tr-tm 

A = -sin{a) e ^2 5m(|3) e ^1 






B^'^ = cos{a)cos{^)e' 

TM TR-TM 



C'«-[(l-e ^1 )co5((3)-l]e ^1 +1 
If we apply the sequence A^^ times and re-arrange the terms we get the series: 

N 
Mf ^^(5) = Mo +Mo [A'^'^cosiqs) +5^^^ + C^^ - 1] £ [A'^^cosiqs) +5^^]"-^ 

(5) 
for the starting magnetization state after A^ repetitions of the sequence. 

Summing an infinite number of terms results in the expression for the steady state 
M^^{s) after a large number of TR periods: 

ee, , A^'^cos(qs)+B'^^ + C^'^-\^ 

One can then calculate the magnetization state after the (3 pulse in the steady state: 

Mf^^{s) = [A^cos{qs) + B^]Mf{s) +C^Mo (7) 

and after TB: 



Mf '^^(5) = [A^'^cos{qs) + B'^^]Mf{s) +C'^'^Mo (8) 

We show graphs of Eq. [6], [7], and [8] in Figure 2 for TR — 2s. 
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4 Results 



We now show in Figure 3 representative Mf^{s) magnitude images obtained with 
the sequence described in section 2 for four different values of TR = 5s, 2s, Is, 500ms. 
Figure 4 shows several cross sections through row #128 of Figure 3. The object 
is an 18mm glass sphere filled with silicone oil. Data points are superimposed 
with the corresponding magnitude of the theoretical curve. The Ti of the silicone 
oil (at 400MHz) was measured by spectroscopic inversion recovery to be 1.4s. A 
Bruker DRX400 Micro 2.5 system was used with a custom 27mm diameter 3 IP/IH 
birdcage coil. 10 TR periods were used to establish steady state. A 10s "relax" 
delay was used between phase encodes to establish full relaxation. Gq was 3ms 
and 2.5mT/mm, with Gspoui of 5ms and lOOmT/mm. No attempt was made to ac- 
count for B\ inhomogeneity. A single scaling parameter was used for all theoretical 
curves. We achieved good agreement with the theoretical predictions. In the se- 
quence as used TM = TB = 1ms. A variety of other Gq directions and strengths 
show similar agreement with theory. Better agreement in the fit between experi- 
ment and theory can be obtained with a = (3 = 75°than with the nominal 90°. A 
B\ map needs to be determined to see if this corresponds more closely to the actual 
experimental conditions. 



5 Conclusions 



The expressions developed and verified above should be useful to those wishing to 
understand or utilize harmonics in the Mf^{s) profile in DDF based sequences in 
the situation where the diffusion distance during TR compared with X in negligible. 
This is especially true for those carrying out structural measurements which depend 
on a well defined correlation distance. The theory should also hold for spatially 
varying magnetization density Mq = Mo{r), and longitudinal relaxation Ti = T\ (r). 
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Figure 4. Row 128 data (points) and 

a = p = 90° , TT? = 2s, TM = TB = 1ms, Ti = 1 As relax = lOs. 
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